Theory of stochastic transitions in area preserving maps 



Massimo Tessarotto 1 a,c and Piero Nicolini 2 a, ~ 



"Department of Mathematics and Informatics, University of Trieste, Italy 
b National Institute of Nuclear Physics (INFN) Trieste Section, Italy 
c Consortium for Magnetofluid Dynamics^ University of Trieste, Italy 

Abstract. A famous aspect of discrete dynamical systems defined by area-preserving maps is the physical interpretation of 
stochastic transitions occurring locally which manifest themselves through the destruction of invariant KAM curves and the 
local or global onset of chaos. Despite numerous previous investigations (see in particular Chirikov [1], Greene [2], Percival 
[3], Escande and Doveil [4] and Mac Kay [5]) based on different approaches, several aspects of the phenomenon still escape 
a complete understanding and a rigorous description. In particular Greene's approach is based on several conjectures, one 
of which is that the stochastic transition leading to the destruction of the last KAM curve in the standard map is due the 
linear de stabilization of the elliptic points belonging to a peculiar family of invariants sets {I(m,n)} (rational iterates) having 
rational winding numbers and associated to the last KAM curve. Purpose of this work is to analyze the nonlinear phenomena 
leading to the stochastic transition in the standard map and their effect on the destabilization of the invariant sets associated 
to the KAM curves, leading, ultimately, to the destruction of the KAM curves themselves. 



INTRODUCTION: PHENOMENOLOGY OF STOCHASTIC TRANSITIONS 

The phenomenology of stochastic transitions for canonical (and area-preserving) maps, such as the standard map [1], 
is one of the most popular subjects of investigation in nonlinear dynamics. In particular, in the case of the transition 
to global chaos this amount to the search of the the so-called last KAM curve, i.e. a continuous invariant curve 
which separates nearby subdomains of the phase-space (F) characterized at most only by local chaos. The search of 
such a KAM curve and the conditions of its possible destruction for arbitrary canonical maps may result potentially a 
difficult problem. In the past, sufficient conditions for the global stochastic transition, for the standard map and suitable 
families of area-preserving maps, have been proposed, such as Chirikov Criterion [1]) or Greene's Chaotic Hypothesis 
[2]. However, despite various attempts (see in particular Chirikov[l], Greene [2], Percival [3], Escande and Doveil [4] 
and MacKay [5]) based on different approaches, several aspects of the phenomenology of stochastic transitions still 
escape a complete understanding and a rigorous description. 

In two dimensions closed KAM curves can be closed (so-called invariant circles) or open. In the first case they 
separate phase space since trajectories cannot cross KAM curves, and hence produce the phenomenon of phase- 
space confinement in phase space. Also open KAM curves, however, can separate space. Discrete invariant sets J 
do not separate space and therefore cannot give rise to phase-space confinement. In general one expects that the 
winding number of an arbitrary closed (or open) KAM curve, in particular those which separate space, depend on the 
stochasticity parameter K, i.e., 

a = a(J,K). 

Therefore, a KAM curve can be destroyed for appropriate values of K, i.e., when a(J,K) becomes rational. This 
is called local stochastic transition of the KAM curve J. The local destruction of a KAM surface, and consequent 
disappearance of local phase-space confinement, is also called transition to local chaos. In particular, there may 
exist a special KAM curve, the so-called last KAM curve, whose destruction permits the existence of phase-space 



1 email: M.Tessarotto@cmfd.univ.trieste.it 
- email: Piero.Nicolini@cmfd.univ.trieste.it 
3 Web site: http://cmfd.univ.trieste.it 



trajectories which "connect" all boundaries lines of phase space T, thus destroying "global" phase-space confinement. 
The stochastic transition of this KAM curve is called as global stochastic transition or transition to global chaos. 



A CLASSIFICATION OF AREA-PRESERVING MAPS 



For the purpose of attempting a systematic investigation of area-preserving maps, generated by discrete dynamical 
systems, a convenient approach is to construct them, in analogy to continuous dynamical systems, in terms of a 
variational principle. Indeed the same concept of Hamiltonian system can be extended to discrete dynamical systems 
in such a way to yield, by construction, area-preserving (i.e., conservative) maps. By analogy to continuous theory, 
discrete Hamiltonian systems can be defined by means of a suitable discrete variational principle expressed by a 
variational equation of the form 8S(x) =0 and denoted as "discrete modified Hamilton variational principle". As 
usual, this equation must be satisfied by arbitrary syncronous variations of the form 5x, = (5q,-, 5p,) (with ;' = 0,n), 
being <5x, = x, — xj, and x,-,xi ( - two arbitrary elements of the syncronous functional class 



{x} = {x,- : x t ,i e N,x = x A ,x„ = x B ,x,- = (q,-,p,) GTCR 2 *} . 



(1) 



where x = x A ,x„ = x B are suitable compatible boundary conditions and q = {q\, ..,q g ),p — (pi,..,p g ), being g an 
integer g > 1 to be identified with the degree of freedom. A convenient definition of the variational functional is 



S(x)= £ {(qi+i-q;)-P;+i-#(q;>Pi+i)}, 

ieO.n 



(2) 



where f/(q,,p ( - + i) denotes a suitably smooth real function, here denotes as discrete Hamiltonian function. The 
corresponding Eulero-Langrange equations read therefore: 



8S{x) 
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(p<-Pi+i)-^— ^(q/,Pi+l) 



= 0, 



where the partial derivatives should be meant in the sense 

d 



H(qi,p i+ i) = —H( qi ,p) 
dpi+l op 

^-//(q ; ,p I+ i) = ^#(q,p i+1 ) 



p=Pi+i 



q=q, 



(4) 



(5) 



(6) 



Hence, they deliver the discrete canonical map T 

' Qi+i = q< + 3j|^#(qj.p<+i). 
p/+i =P/-^^(*''P«'+ 1 )' 

X —X , X n = X , 

(defined for ;' = 0,n) which can be proven to satisfy the fundamental Poisson brackets: 

[q<+l.P<+l](q i ,p l )=l. ( ? ) 

[q;+i,qi+i]( q „ P ,.) = [Pi+i,Pi+i]( q ,, Pi ) = »• (8) 

This implies, by virtue of the previous definition (2) for the variational functional, that discrete Hamiltonian systems, 
just like continuous ones, result trivially area-preserving. In particular, this property is satisfied by arbitrary elements 
of the sequence {xj, i € N} x a , defining the phase-space trajectory of the discrete canonical system, which is determined 
by the map T , for an arbitrary initial state x = x A e T (and letting x„ = x B ). 

Conversely, in the particular case of two-dimensional (i.e., forg = 1), it follows that area-preserving maps T of the 
form 

*<= (ihPi ) -» *i+i={<li+vPi+\) = T *i, (9) 
are necessarily canonical and hence of the type defined above (6). 



Parameter-dependent canonical maps 



In the sequel we shall consider, without loss of generality, parameter-dependent 2D canonical maps T charac- 
terized by a bounded phase space T = [— n,n[x [0,2n[ and having a discrete Hamiltonian function of the form 
H = H(qi,pi + i,K), being K a real parameter (to be denoted as stochasticity parameter) and H a real analytic function 
of K. Moreover, we shall impose that the discrete Hamiltonian H satisfies the constraint 

H( qi ,p i+1 ,K = 0)=H (p i+ i), (10) 
for arbitrary i = 0,n and phase-space trajectories, which implies 

CO 

H{q ilPi+u K) = H ( Pi+l )+ £ KiHj( qi ,p i+l ). (11) 
An example of 2D maps of this type, which results also a twist map, i.e., satisfies the twist condition 



dpi+1 

is provided by the so-called standard map [[1],[2]]. 



t^O, (12) 



= qi+Pi+i (13) 

Pi+l = pi+Ksinqi. 

Since the map is manifestly area-preserving, it results also canonical, the discrete Hamiltonian function being in this 
case: 

H(q h p i+1 ,K) = -pj +1 +Kcos qi . (14) 

Two-dimensional parameter-dependent area-preserving maps of this type are characterized by discrete and continuos 
invariant sets (the latter called KAM curves). In particular, for K — all continuous subsets of T defined by the 
equations 

Pi=p A (15) 

are necessarily invariant (KAM curves). For K > KAM curves can change shape or disappear altogether. 

In particular, discrete invariant sets may appear which are characterized generally by a finite number (n denoted as 
order of the set) of m-periodic invariant points (i.e., points which are returned to themselves after m iterations of the 
map). Discrete and continuous invariant sets can be distinguished by means of their winding numbers a{x A ), i.e.: 

fiW (k a ) 

a{x A ) = lim - — i — '- (16) 
where 5^ (a(x A ) is the rotation angle after at the k— th iteration of the map T, namely: 

s ( V)= E (pj^-pj) d7) 

7=1,*- 1 

and xj is the j— th image of x A through the canonical map T. As a consequence, it is immediate to prove that discrete 
invariants sets of order n and periodicity m, to be denote by the symbol I(m,n) = { m,n,K,x A } , display necessarily 
a rational winding number a(x A ) = m/n, whereas KAM curves have necessarily irrational winding numbers. In both 
cases, one expects generally that winding numbers result function of the stochastic parameter K, i.e. that a = a(J,K), 
where J denotes an arbitrary invariant set. In particular one may expect that a given KAM curve y which exists for 
K = 0, and is characterized by an irrational winding number, may become a discrete invariant set in correspondence of 
a suitable critical values of the stochasticity parameter K the function a(J,K) may become rational. This is called 
local stochastic transition. A basic conjecture (Greene, 1979 [2]) is that for arbitrary area-preserving maps each 
KAM curve y can actually be constructed as a limit set of an appropriate family of m-periodic discrete invariant 
sets {I(m,n),m = m(n),n e , i.e., forn — > oo ; 

I(m(n),n)^Y, (18) 




FIGURE 1. Location of last KAM curve and of the rational iterates 



where the discrete invariant sets I(m(n),n) are characterized by the distinctive property that their winding numbers 
a(I(m,n)) — tend to a(y) in the asymptotic limit, i.e., 

m(fi) 

Urn a{I{m,n))= lim — ^ = a(y). (19) 

A basic conjecture proposed by Greene [2]) is that for arbitrary area-preserving maps each KAM curve y, having 
winding number a(y), can actually be constructed as a limit set of an appropriate family of m-periodic discrete 
invariant sets (to be denoted as rational iterates) {I(m,n),m = m(n),n £ N} , i.e., for n — > °° : 

7(m(»),«)^y, (20) 

where the discrete invariant sets J(m(n),n) are characterized by the distinctive property that their winding numbers 
a{I{m,n)) — tend to a(y) in the asymptotic limit, i.e., 

m(n) 

lim a(I(m,n)) = lim — — = a(y). (21) 



GREENE'S CHAOTIC HYPOTHESIS 

The method developed by Greene [2] for investigating stochastic transitions in area-preserving maps (and in particular 
for the standard map) is based on the following set of conjectures, denoted as Chaotic Hypothesis: 

A) 1st conjecture (irrational winding number of the last KAM curve): the last KAM surface is characterized by the 
irrational winding number 

Vs-i 

«(y) = «gm = — 2 — ; {2T > 

{Uqm = golden mean) which can be represented by the continuous fraction 



B) 2nd conjecture (mechanism of destruction of a KAM curve) : the disappearance of a KAM surface (in particular 
the "last KAM surface" whose destruction produces the global stochastic transition) is related to the linear destabi- 
lization of the invariant points (elliptic equilibrium points) of an appropriate family of m-periodic discrete invariant 



sets I(m,n), i.e. when the Lyapunov exponents of the invariant points become positive for suitable values of the 
stochasticity parameter K . The stochastic transition can be estimated by analyzing the bifurcation conditions (linear 
destabilization) of the invariant sets I(m,n) for sufficiently high values of the order n of these sets. In particular, for 
sufficiently high values of n and for all n \ > n all points of the sets I(m{n\),n{) become linearly unstable. 

C) 3rd conjecture (invariant sets): in the case of the last KAM surface of the standard map the discrete invariant sets 
J(m, n) coincide with the so called rational iterates. The winding numbers of these special invariant sets are given by 
the following sequence, obtained truncating the previous continuous fraction, i.e., 

i 1 in(ii) n 

a(/(l,2)) = — — = -,...,a(/(m,n)) = -^,a(/(n,n+m)) = -> Ogm, (24) 

1+12 n n + m 

yielding the family of sets 

{/(m(n),n),neiV} = {/(l,2),/(2,3),/(3,5),...,/(in(n),ii),...} (25) 

where n,m are Fibonacci numbers belonging to the sequence 

n,/ne{2,3,5,8,13,...Jt}, (26) 

with k given by the sum of the previous two elements of the sequence. Here m and n denote, respectively, the 
periodicity and the order (number of points) of the generic set of the family I(m,n). 

Using numerical experiments Greene was able to justify at least in part there conjectures and achieve the following 
results for the standard map: 

I) identification of the last KAM curve and the family of rational iterates: (see figure 2){I(m,n),n EN} (see 
figure 5): near stochastic transition the invariant sets I(m,n) are found to be localized below the last KAM; these sets 
continue to exist also after transition (i.e., for values of the stochasticity parameter K larger than the critical value 

Kcrit )• 

II) linear stability analysis of rational iterates I(m,n) : highest order n considered for the rational iterate n = 233, 
i.e., I(m,n) =7(144,233) 

III) estimate of the critical value of the stochasticity parameter K at stochastic transition: 

K = 0.9716, (27) 

determined by the requirement that the elliptic invariant points of the set 7(144,233) become linearly unstable. 

Although these results have been confirmed by other authors using different approaches (Percival[3], Escande and 
Doveil[4] and MacKay[5]), they leave unexplained several aspects and basic issues of the theory. In particular, the 
identification of the invariant sets I(m,n) (conjecture B) associated to the last KAM curve remains not obvious even 
for the standard map. In the same way this problem remains unanswered for arbitrary KAM curves and more general 
area-preserving maps. 



THE NONLINEAR CHAOTIC HYPOTHESIS (NCH) 

In this regard, a natural question is whether the family of invariant sets I(m,n) (the rational iterates), associated to a 
given KAM curve, i.e., for which the KAM curve is a limit set, is unique or not. This issue arises immediately when 
trying to relate Greene's chaotic hypothesis to the well known Chirikov criterion for the overlapping of resonances 
(Chirikov Criterion [1]). In fact, the latter states that KAM curves are destroyed when the distance between the two 
lowest order resonances which encompass the same KAM curve, i.e., lie on opposite sides of the KAM curve, becomes 
comparable to the sum of the semi-amplitudes of the separatrix islands formed by the two resonances. On the other 
hand the existence of multiple sequences of discrete invariants sets having the same limit set (a given KAM curve) 
does not contradict the original Greene conjecture, which does not specify the precise location of the invariant sets 
J(m,n) nor exclude their possible non-uniqueness. This has motivated us to search for possible candidates for these 
invariant sets, associated in principle to an arbitrary KAM curve, which are denoted here as alternate rational iterates 
Y(mj,n). Numerical experiments performed for the stochastic transition of the last KAM curve of the standard map 
[6] confirm this conjecture and suggest a new formulation of the Greene's Chaotic Hypothesis, replacing conjectures 
B and C with the following three new conjectures Bl and CI: 

Bl) new 2nd conjecture (invariant sets related to a KAM curve) : the disappearance of a KAM surface y (in particular 



the "last KAM surface" which determines with its destruction the global stochastic transition) is related to the existence 
to two families of discrete invariant sets I(m,n) and Y(mj,n) , with j a suitable integer, defined in such a way their 
winding numbers tend to a (7), i.e., 

nidi) 

lim a(I(m,n))= lim — ^ = a(y), (28) 

lim a(Y(m, ,«)) = lim m ^ = a(y). (29) 

CI) new 3rd conjecture (invariant sets related to the last KAM curve of the standard map) : in the case of the last 
KAM surface of the standard map the discrete invariant sets I(m,n) and Y(rrij,n) , coincide, respectively, with the 
sets of rational and alternate rational iterates when letting 7 = 1. 

Dl) new 4th conjecture (nonlinear mechanism of destruction of a KAM curve): the stochastic transition of the 
last KAM surface occurs when the Euclidean distance between corresponding elliptic points of the sets I{m,n) and 
Y(mj,n) 7 considered as a function of the stochasticity K, becomes minima, for sufficiently high values of n. 

The set of hypotheses A (Greene's first conjecture), together with Bl.Cl and Dl, is here denoted here as Nonlinear 
Chaotic Hypothesis (NCH). We stress that the conjecture Dl (new 4th conjecture) appears consistent with (but differs 
from) Chirikov Criterion, since the distance between nearby resonances cannot exceed the sum of the semi-amplitudes 
of the separatrix islands formed by the two resonances. Moreover, it replaces the assumption on linear stability with a 
non-linear destabilization criterion, i.e. the requirement that or the Euclidean distance between adjacent elliptic points 
belonging respectively to two sets I(m, n) and Y(m\,n) having the same order n. In fact, in principle it is not obvious 
at all why a linear stability criterion (as Greene's 4th conjecture) should hold in general since an elliptic invariant 
point might result at the same time linearly unstably and nonlinearly stable. In addition, invariant points of different 
invariants sets (for example adjacent points invariants sets of the same order n, I(m,n) andF(mi,n)), may not become 
all linearly unstable simultaneously, i.e. for the same value of the stochasticity parameter K. 

Therefore, NCH should permit, in particular, not only a more accurate evaluation of the critical value of the stochastic- 
ity parameter K occurring at stochastic transition, but also a precise identification of the families of discrete invariant 
sets associated, in principle, to an arbitrary KAM curve. Finally, the new approach closes the gap between Greene's 
Chaotic Hypothesis and Chirikov criterion of overlapping of resonances, permitting at the same time a deeper under- 
standing of the mechanism of stochastic transitions in area-preserving maps. 
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